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(J ; ABSTRACT 

^' 

I ■ Despite decades of study, it remains unclear whether there are distinct radio-loud 

i-( ■ and radio-quiet populations of quasi-stellar objects (QSOs). Early studies were limited 

^ ■ by inhomogeneous QSO samples, inadequate sensitivity to probe the radio-quiet popu- 

lation, and degeneracy between redshift and luminosity for flux-density-limited samples. 
Our new 6 GHz EVLA observations allow us for the first time to obtain nearly complete 
(97%) radio detections in a volume-limited color-selected sample of 179 QSOs more lu- 
VO . minous than Mj = —23 from the Sloan Digital Sky Survey (SDSS) Data Release Seven in 

^r^ . the narrow redshift range 0.2 < z < 0.3. The dramatic improvement in radio continuum 

[-«^ . sensitivity made possible with the new EVLA allows us, in 35 minutes of integration, to 

2 ■ detect sources as faint as 20 /iJy, or log[L6GHz(WHz~^)] ^ 21.5 at z = 0.25, well below 

the radio luminosity, log[L6(WHz~^)] fn 22.5, that separates star-forming galaxies from 
radio-loud active galactic nuclei (AGNs) driven by accretion onto a super-massive black 
^^ . hole. We calculate the radio luminosity function (RLE) for these QSOs using three con- 

straints: (a) EVLA 6 GHz observations for log[L6(WHz-i)] < 23.5, (b) NRAO-VLA 
Sky Survey (NVSS) observations for log[Le(WHz^^)] > 23.5, and (c) the total number 
of SDSS QSOs in our volume-limited sample. We show that the RLE can be explained 
as a superposition of two populations, dominated by AGNs at the bright end and star 
formation in the QSO host galaxies at the faint end. 
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INTRODUCTION 



Since the discovery of optically selected QSOs (ISandagd Il965l ). the difference between the 
"radio-loud" (RLQ) and "radio-quiet" (RQQ) QSO populations remains elusive. There has been 
a continuing controversy as to whether the rad io luminosity distribution of QSOs i s bimodal (e.g., 



Kellermann et al 



200^ 



3|; iLacy et al. 



19891 : llvezic et al.ll2002l . 120041 ) or is merely broad and smooth (e.g.. lCirasuolo et al. 



20011 ). A bimodal distribution suggests that two distinct physical processes are 



present, with one process being significantly more powerful than the other. 

The two primary sources of radio emission from galaxies are (1) accretion onto super-massive 
black holes in active galactic nuclei (AGNs) and (2) star formation. Radio-AGN emission is 
due to synchrotron from relativistic plasma jets, and the associated hotspots and lobes result- 
ing from jet interaction with the surrounding medium. Star formation yields free-free emission 
from HII regions and synchrotron radiation from relativistic electrons believed to be accelerated 
in supernova remnants. Radio AGNs can be ex tremely luminous , with radio emission reaching 
log[L6GHz(WHz~^)] ^ 27 in the case of 3C 273 f lKellermann et al.lll969l ). whereas star formation 
leads to radio luminosities of log[L6GHz(WHz~^)] fa 21 for a galaxy like the Milky Way, or as high 
as log[L6GHz(WHz~^)] ^ 23.5 in the case of a star-bursting galaxy such as Arp 220. Radio- AGNs 
dominate the radio luminosity function (RLF) of galaxies brighter than log[L6GHz(WHz'^)] ^ 22.5 
i n the radio, while star-formation in galaxies without an AGN dominates at fainter luminosities 
JCondon et al.ll2002h . 



The question as to the nature of radio emission in radio-quiet QSOs has not been previously 
addressed using a homogeneous and complete QSO sample that is sensitive to luminosities sig- 
nificantly fainter than log[Lj.adio(WHz~^)] = 22.5. It is possible that some, if not all, QSOs are 
hosted in star-forming galaxies, which only become visible at radio wavelengths when the AGN 
radio emission is below sor ne given threshold . In this Letter, we use results from the Expanded 
Very Large Array (EVLA) (jPerley et al.ll201ll ) to constrain the radio luminosity function (RLF) of 
QSOs. We show that the RLF is consistent with the hypothesis that QSO radio sources with 6 
GHz spectral luminosity log[Le(WHz~^)] > 23 are powered primarily by AGNs, while those with 
log[L6(WHz~^)] < 23 are powered primarily by star formation in their host galaxies. 

By taking advantage of recent advances in both radio and optical capabilities, we have, for the 
first time, obtained nearly complete radio detections in a large volume-limited sample of optically 
selected QSOs. In this Letter, we describe early results of our EVLA observations. In Section [2] 
we describe the target selection and radio observations. In Section [3] we present our results and 
interpret them in terms of the superimposed radio contributions from AGNs and star-forming host 
galaxies. Our conclusions are summarized in Section HI 
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DATA 



2.1. Target selection 



Our color-selected targets were drawn from the ISchneider et al.l ( 20101) quasar catalog of the 
seventh data release (DR7) of the Sloan Digital Sky Survey (SDSS) (JAbazajian et al.ll2009l ). We 
selected nearby QSOs in order to reach a spectral luminosity limit of log[L6(WHz~^)] ^ 21.5 after 
35 minutes of integration with the 2 GHz bandwidth available at the EVLA, and in a narrow redshift 
range to minimize the effects of evolution. Our selection criteria were: 



1. targeted by the SDSS for spectroscopy using the "low z" color criteria (see explanation below) 

2. Mi < -23 {Ho = 70 km s'^ Mpc-\ ^a = 0.7) 



3. 14 < z < 19 (apparent magnitude i corrected for extinction according to lSchlegel et al.lll998l ). 
which is well above the SDSS completeness limit 

4. 0.2< z < 0.3 

5. 6 > 30° 



We used the SDSS to select a homogeneous sample of SDSS QSOs. The SDSS identifies low- 
redshift QSO candidates for follow-up spectroscopy by their non-stellar colors in the SDSS ugri 
color cube. These sources are labeled in the SDSS database with the QSO targeting flag "low z". 
QSOs are then selected from all objects with spectra, and identified by ha ving at least one broad 
emission line with FWHM greater than 1000 km s~^ (JRichards et al.ll2002l ). 



The selection criteria result in 179 targets that comprise a volume-limited, color-selected sample 
of all Mi < —23 QSOs in ~2.66 steradians in the redshift range 0.2 < z < 0.3. Owing to the color 
selection criteria, our sample may miss some QSOs in our search volume that happen to have stellar 
colors. However, it constitutes a fully optically selected sample with no radio biases. About one- 
third of QSOs in our target list were already known rad io emitters found in the Faint Images of the 
Radio Sky at Twenty cm (FIRST) (JBecker et al.l 119951 ) survey (peak flux density Sp > 1 mJy at 
Y = 1.4 GHz) and/or found in the NRAO— VLA Sky Survey (NVSS) {Sp > 2.5 mJy at z/ = 1.4 GHz) 
JCondon et al.lll998l ). 



2.2. Observations 

All of the observations were made with the C configuration of the EVLA using a nominal 
bandwidth of 2 GHz centered on 6 GHz in each of two circular polarizations, which were combined 
to make total intensity images. Because any source detected in the FIRST or NVSS surveys is 
expected to be at least a few hundred microjanskys at 6 GHz, we observed the QSOs identified 
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with FIRST or NVSS radio sources for only three minutes. To minimize sidelobes in the dirty beam 
we tried to divide these short snapshot observations into two 90 s segments widely separated in hour 
angle, but owing to scheduling constraints about 10% of these FIRST/NVSS sources were observed 
for only 90 s at one hour angle. All such sources were detected in the single 90 s observation. All 
other QSOs were, at first, observed for five minutes. Most of those that were clearly detected in 
the five minute observation were re-observed for another five minutes at a different hour angle. The 
remaining sources were observed for an additional 30-35 minutes, their observations were combined, 
and the resulting images generally reach the expected rms noise levels of 6-8/iJy per synthesized 
beam solid angle (~ 3'.'5 FWHM). 

For the purposes of this Letter, we consider sources to be detected at a 3a confidence level. 
This confidence level is appropri ate, given that the QSO positions are known with ~ O'.'l accuracy 
from the SDSS fjPier et al.ll2003[ ). Less than 10% of the QSOs in our sample have 3-4cr detections; 



the probability than any one of these is false is < 0.0014. Thus the probability that our sample 
contains even one false detection is < 2.5%. 



3. AGN and Galaxy Components of QSO Radio Emission 

In the top panel of Figure 1, we show the distribution of observed 6 GHz flux densities for the 
compact (unresolved on ~ 3'.'5 scale) radio components coincident with the 179 QSOs in our color- 
selected sample. In the bottom panel, we show the corresponding spectral luminosity distribution 
calculated assuming a spectr al index a = —0.7 (where L^, oc u""; see discussion below) for the K- 



correction ( iHogg et al.ll2002l ). The observed distributions are clearly peaked around 100 /iJy and 
a few times 10^^ WHz~^ respectively, which are well above our detection limits. Only six QSOs 
remain undetected. For one of these, 103421.71+605318.1, there was a strong 69 mJy source in the 
field, which degraded the sensitivity limit by a factor of 2.5. For the other five undetected QSOs, 
the rms noise ranges from 7 to 11 fiJj, somewhat above our nominal limit owing to interference, 
limited on-source observing time, or inadequate dynamic range. We include the 3a upper limits of 
these six undetected sources in Figure 1. 

For the purposes of this Letter, we assume that all of the QSOs have a spectral index of 
a = —0.7 between 6 GHz and 1.4 GHz. A direct measurement of the true spectral index for a faint 
FIRST source is difficult, as a signific ant fraction of emission can be missed from faint sources owing 
to "clean bias" ( Condon et al.lll998l ). The NVSS resolution (~ 45" beam) is not a good match to 



the resolution of our EVLA observations (~ 3'.'5 beam); thus it can also be difficult to measure 
accurate spectral indices for extended NVSS sources in the sample. The 6 GHz EVLA observations 
suggest that, except for the very strongest QSOs, the sources appear unresolved. Therefore, we can 
use NVSS images to estimate the 1.4 GHz flux density of the remaining QSOs, even though they do 
not have counterparts in the NVSS catalog. We measured the flux densities at the optical positions 
of all QSOs not identified with cataloged NVSS sources stronger than the catalog limit of 2.4 mJy. 
The fiux-density distribution of these QSOs has a median value of 330 +- 30 microJy. The median 
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Fig. 1. — Top: Histogram of the distribution of 6 GHz flux densities observed by the EVLA in our 
sample of 179 QSOs. For six undetected sources, 3a upper limits are shown (arrows). Bottom: 
Intrinsic spectral radio power of the 179 QSOs at 6 GHz in the source frame. 
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flux density of our 6 GHz EVLA detections is 126 /iJy, yielding an estimate of a = —0.6 9 for the 



avera ge spectral index of the QSO sample. This value is typical of star-forming galaxies (ICondon 



19921 ). Using the actual spectral indices of individual sources would not significantly change Figure 1 
nor the analysis discussed in the remainder of this section, as a change in spectral index of Aa = 0.5 
corresponds to a factor of two in luminosity and a difference of only 0.3 on a log scale. 

The 6 GHz spectral luminosity function of 0.2 < z < 0.3 color-selected QSOs derived from 
these observations is shown by the black data points with error bars in Figure 2. T he luminosity 



function was calculated using the \^/\4iax rnethod on detected sources (jSchmidtlll968l ). The V"/l^ax 
method statistically determines the RLF using detected quantities; upper hmits are not used in the 
calculation. Instead, undetected sources are accounted for by the normalization of the accessible 
volume iy) of the detected sources. Because there are only six undetected QSOs, there cannot 
be more than six undetected QSOs in the luminosity range 18.8 < log[L6(WHz^^)] < 21.2, which 
fact we indicate by the wide upper- limit symbol in Figure 2. Also shown are the 6 GHz luminosity 
functions of nearby {z < 0.05) galaxies whose radio luminosities are dominated either by AGNs (solid 
green curve) or star formation (solid red curve), calculated from the 1.4 GHz luminosity functions 



(jCondon et al.l |2002| ) using spectral index a = —0.7. The 6 GHz spectral luminosity function of 
nearby galaxies whose radio sources are dominated by AGNs extends to higher radio luminosities 
than the luminosity of star-forming galaxies, so AGN-powered sources are more common above 
log[L6(WHz~^)] i^ 22.5; star-forming galaxies have higher space densities than AGNs at lower 
luminosities. The nearby "starburst" galaxy M82, which is a typical radio-selected star-forming 
galaxy, has log[Lfi(WHz~ ^)] ^ 21.5. This va lue is close to the 10^^ W Hz~^ spectral luminosity of 



the Milky Way at 6 GHz flBerkhuiisenl 119841 ) 



Above log[L6(WHz'^)] ^ 23.5, the spectral luminosity of the ultraluminous starburst galaxy 
Arp 220, it is likely that the QSO radio emission is primarily produced by AGNs. To estimate the 
luminosity function that would result if all QSO radio emission were powered entirely by AGNs, 
even at lower luminosities, we extrapolated the high-luminosity QSO luminosity function to faint 
luminosities as shown by the dashed green curve in Figure 2 for 21 < log[L6(WHz^^)] < 25.5 using 
a calculated slope of —0.30 ±0.03. The slope for this extrapolation is not adequately constrained by 
the small number of log[L6(WHz~^)] > 24 QSOs in our sample (16 of 179). In stead, we determined 



the slope from the observed 1.4 GHz QSO luminosity function, using NVSS (ICondon et al.lll998[ ) 
measurements at the positions of SDSS QSOs in the redshift range 0.2 < z < 0.45. The lack of 
observed QSOs with very high radio luminosities motivates the fall-off at log[L6(WHz^^)] > 26, 
although its exact form is not known and cannot be determined from our data. However, this fall-off 
does not affect the calculation of the power-law slope below log[L6(WHz^^)] ^ 25.5. The slope 
is consistent with the EVLA data points for log[L6(WHz~^)] > 24, as it should be. The dashed 
green curve is shown extrapolated to log[L6((WHz)~ )] ^ 21.3 to emphasize that it falls below the 
observed luminosity function in this range; either the AGN luminosity function has a "bump" not 
observed in the AGN luminosity function of nearby galaxies or there is an additional energy source 
contributing to the radio emission of most QSOs with 21 < log[L6(WHz~^)] < 23. If the luminosity 
function of the AGN component alone is extrapolated to still lower luminosities (as suggested by 
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the dotted green line), it must fall off around log[L6(WHz^^)] ~ 19 lest the number of AGN exceed 
the total number of QSOs in our SDSS sample. The exact form of the luminosity function and its 
cutoff indicated by the dotted curve is unknown but not critical, because it always predicts that 
most QSOs powered only by AGN would be weaker than the EVLA detection limit S* ~ 20 /xJy, or 
log[L6(WHz-i)] > 21.3 if ^ > 0.2. 

In order to be consistent with both the EVLA detections of QSOs stronger than 20 /iJy and the 
total number of SDSS QSOs, the QSO RLF must rise sharply just below log[L6(WHz~^)] ^ 23.5 
and fall fairly sharply at lower luminosities. The total number of SDSS QSOs in this luminosity 
function is known (179 for 0.2 < z < 0.3j, and limits the integral area of the RLF. Even if all 
six non-detections are not in the luminosity range 18.8 < log[L6(WHz~^)] < 21.2, the upper limit 
in Figure 2 constrains the most conservative estimate of the RLF fall-off. If any of the six non- 
detections actually have radio luminosity in the range 21.2 < log[L6(WHz^^)] < 22, the RLF at 
low luminosities must fall even more sharply than what is shown in Figure 2. 

We suggest that the rise in the RLF at log[L6(WHz^^)] ^ 23.5 is the result of radio emission 
from the star-forming host galaxies. The solid red curve in Figure 2 show s the 6 GHz space d ensity 



of radio sources powered by all star- forming galaxies in the local universe (jCondon et al.ll2002l ). The 
dashed red curve illustrates one possible form of the 6 GHz spectral-luminosity function for the host 
galaxies of 0.2 < z < 0.3 QSOs, a parabolic fit with fall-off determined independently on either side 
of the peak. The peak of the curve is at log[L6(W Hz~^)] = 22.2 and log[pm(Mpc~^ mag~^)] = —7.2. 
The slope of the fall-off on either side of the peak is not tightly constrained by our data. The fall-off 
slope indicated by of our dashed red curve at log[L6(WHz~^)] > 22.5 is consistent with the fall-off 
slope of the solid red curve, but with 3% of the space density. This result suggests that the space 
density of luminous starbursts in QSOs at 0.2 < z < 0.3 is approximately 3% the space density 
of local galaxies powered by starbursts. However, this fraction may be different if the true RLF 
distribution of star-forming QSO host galaxies is not well-represented by the dashed red curve. 

The black curve is the corresponding luminosity function of QSOs whose radio luminosities 
are the sum of both the AGN and star-forming galaxy luminosities, on the assumption that the 
AGN and star-forming luminosities are statistically independent. The spectral-luminosity functions 
indicated by the dashed red curve and by the black curve are constrained by the EVLA detections 
shown as black data points in Figure 2 and by the total number of SDSS QSOs. The peak of the 
QSO fiux-density distribution implied by the black curve is (log[S'(/iJy)]) ^ 2.1 at z/ = 6 GHz, and 
the peak of the the QSO host-galaxy luminosity function should be at log[L6(WHz^^)] ^ 22.4. 

As it must, the RLF of the QSO host galaxies (dashed red curve) lies below the RLF of all 
galaxies powered primarily by star formation (solid red curve). The dashed red curve lies at the 
luminous end of the star-forming galaxies RLF, suggesting that most QSO host galaxies have higher 
star-formation rates than galaxies without a QSO. This result is surprising, given that QSOs are 
traditionally associate d with "red and dead" elliptical gala xies, at least for radio-luminous QSOs 



(typically "quasars") (JDunlop et al.ll2003l : iFloyd et al.ll2004l ). If the host galaxies of all QSOs were 



similar to the massive ellipticals with low star-formation rates that typically host radio-luminous 
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Fig. 2. — Our models for the spectral luminosity func tions of 0.2 < z < .3 QSOs are compared 
to the 1.4 GHz luminosity function of nearby galaxies (jCondon et al.l 120021 ) transformed to 6 GHz 
with an assumed spectral index of -0.7. Radio sources powered primarily by star formation are 
shown with a solid red curve and those powered by AGNs as a solid green curve. Solid black points 
correspond to our EVLA data; open blue points correspond to NVSS sources. The dashed green 
curve is an extrapolation of the high-radio-luminosity QSOs to low luminosities, using the slope 
determined from the NVSS data. The dashed red curve represents the spectral luminosity function 
of QSO hosts that are powered primarily from star formation, constrained by the EVLA data. The 
upper limit symbol at log[L6(WHz"^)] ~ 20 represents the six EVLA non-detections. While it is 
not certain that all of the non-detections have luminosities this low (see Figure 1), this upper limit 
constrains the most conservative estimate of the RLE fall-off (see text for details). The black curve 
is the calculated luminosity function for QSOs whose radio sources are powered by both AGNs and 
star formation in their host galaxies. The total area under the black curve is constrained by the 
number of SDSS QSOs in the volume-limited (0.2 < z < 0.3) sample. 
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QSOs, we would expect the RLF of our QSO sample to follow the extrapolation of the dashed 
green curve in Figure 2. The rise in the RLF coincides with the typical level of radio emission from 
star-forming galaxies, suggesting that the host galaxies of the log[L6(WHz~^)] ~ 22.4 QSOs are 
not "red and dead" elliptical galaxies. Instead, these host galaxies may have spiral morphology, or 
the star formation may be the result of galaxy mergers. It has been suggested that active galaxies 
populating the "green valley" , with colors intermediate between the more typical blue or red galaxy 
colors, may be high-Eddington-ratio ellipticals with some star formation, or the products of mergers 



between mass i ve sp heroidal galaxies and less massive gas-rich galaxies (e.g., ISchawinski et al.ll2010l : 



Kaviraj et al.l 120091 ) . Characterizing the host galaxy properties of these QSOs will be a crucial 



follow-up step in investigating this result. 



4. CONCLUSIONS 

For the first time, our observations adequately sample the radio-quiet population; that is, they 
detect nearly all optically selected QSOs in a volume-limited sample by reaching log[L6(WHz~^)] = 
21.5. Earlier studies did not have the sensitivity needed to study the full range of radio-quiet QSOs 
by reaching faint radio flux densities (S* ^ 1 mJy) for optical samples with significant numbers of 
low-redshift {z < 0.5) QSOs to probe the low end of the RLF. Analyses based on radio-selected 
samples are inherently biased toward the radio-loud population, but previous studies reporting 
smooth luminosity distributions did not extend sufficiently faint to study the radio-quiet population. 

The 6 GHz RLF of low-redshift color-selected QSOs is constrained by EVLA detections of 
sources stronger than 20/iJy and by the very small fraction (6/179) of nondetections in the SDSS 
sample. The strong sources constrain the AGN contribution above log[L6(WHz~^)] ^ 23.5. The 
fainter sources imply a second radio contributor narrowly peaked aroun d log[Lfi(WHz~^)] ^ 22. 4 



and confirm the two-population model of QSO fiux densities found by iKellermann et al.l (119891 ). 
Using the RLFs of nearby galaxies as examples, we suggest that this second contributor is radio 
emission produced by strong star formation in the QSO host galaxies, most of which therefore 
cannot be "red and dead" massive ellipticals. 

In this Letter, we have presented the preliminary results of our investigation to characterize 
the radio population of QSOs using an optically selected, volume limited sample. We suggest that 
the RLF is a superposition of radio emission from AGNs and from host star formation. It will be 
important to better understand the relation between QSO emission in different wavelength regimes 
by characterizing the host galaxies of these QSOs using optical and infrared colors and morphology, 
as well as the radio morphology. The full details of this study will be presented in a separate paper. 
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proved the clarity of the paper. The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by Associated Universities, 



-10- 

Inc. Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, 
the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, 
the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max 
Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site 
is |http://www.sdss.org/| 

REFERENCES 

Abazajian, K. N. et al. 2009, ApJS, 182, 543 

Becker, R. H., White, R. L., & Helfand, D. J. 1995, ApJ, 450, 559 

Berkhuijsen, E. M. 1984, A&A, 140, 431 

Cirasuolo, M., Celotti, A., Magliocchetti, M., & Danese, L. 2003, MNRAS, 346, 447 

Condon, J. J. 1984, ApJ, 287, 461 

— . 1992, ARA&A, 30, 575 

Condon, J. J., Cotton, W. D., & Broderick, J. J. 2002, AJ, 124, 675 

Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, Q. F., Perley, R. A., Taylor, G. B., & Broderick, 
J. J. 1998, AJ, 115, 1693 

Dunlop, J. S., McLure, R. J., Kukula, M. J., Baum, S. A., O'Dea, C. P., & Hughes, D. H. 2003, 
MNRAS, 340, 1095 

Floyd, D. J. E., Kukula, M. J., Dunlop, J. S., McLure, R. J., Miller, L., Percival, W. J., Baum, 
S. A., & O'Dea, C. P. 2004, MNRAS, 355, 196 

Hamilton, T. S. 2010, MNRAS, 407, 2393 

Hogg, D. W., Baldry, I. K., Blanton, M. R., & Eisenstein, D. J. 2002, ArXiv Astrophysics e-prints 

Ivezic, Z. et al. 2002, AJ, 124, 2364 

Ivezic, Z. et al. 2004, in Astronomical Society of the Pacific Conference Series, Vol. 311, AGN 
Physics with the Sloan Digital Sky Survey, ed. G. T. Richards & P. B. Hall, 347 

Kaviraj, S., Peirani, S., Khochfar, S., Silk, J., & Kay, S. 2009, MNRAS, 394, 1713 

Kellermann, K. I., Pauliny-Toth, I. I. K., & Williams, P. J. S. 1969, ApJ, 157, 1 

Kellermann, K. I., Sramek, R., Schmidt, M., Shaffer, D. B., & Green, R. 1989, AJ, 98, 1195 



- 11 - 

Lacy, M., Laurent-Muehleisen, S. A., Ridgway, S. E., Becker, R. H., & White, R. L. 2001, ApJ, 
551, L17 

Pier, J. R., Munn, J. A., Hindsley, R. B., Hennessy, G. S., Kent, S. M., Lupton, R. H., & Ivezic, 
Z. 2003, AJ, 125, 1559 

Perley, R. A., Chandler, C. J., Butler, B. J., & Wrobel, J. M. 2011, ApJ, in press 

Richards, G. T. et al. 2002, AJ, 123, 2945 

Sandage, A. 1965, ApJ, 141, 1560 

Schawinski, K. et al. 2010, ApJ, 711, 284 

Schlegel, D. J., Finkbeiner, D. R, & Davis, M. 1998, ApJ, 500, 525 

Schmidt, M. 1968, ApJ, 151, 393 

Schneider, D. R et al. 2010, AJ, 139, 2360 

Stawarz, L. 2010, in Astronomical Society of the Pacific Conference Series, Vol. 427, Accretion 
and Ejection in AGN: a Global View, ed. L. Maraschi, G. Ghisellini, R. Delia Ceca, & 
F. Tavecchio, 357 

Stawarz, L., Sikora, M., & Lasota, J. 2008, in Astronomical Society of the Pacific Conference 
Series, Vol. 386, Extragalactic Jets: Theory and Observation from Radio to Gamma Ray, 
ed. T. A. Rector & D. S. De Young, 169 



This preprint was prepared with the A AS lA'QTpC macros v5.2. 



